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We study the single-plane-wave Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) states for BCS-Bose-Einstein 
condensation (BEC) crossover at general temperatures T. Because we include the important effects of non- 
condensed pairs, our T 7^ phase diagrams are different from those reported in earlier work. We find that 
generalized LOFF phases may be the ground state for a wide range of (weak through moderately strong) inter- 
actions, including the unitary regime. However, these LOFF phases are readily destroyed by non-zero T. 



Recent atomic physics discoveries in the field of ultracold 
polarized yj, y, |3| fermionic superfluids have important im- 
plications in color superconducting quark matter as well as in 
dense nuclear matter Holla]- Moreover, there has been a long 
standing interest from the condensed matter community Ij in 
observing the very elusive Larkin-Ovchinnikov-Fulde-Ferrell 
(LOFF) [8] states of a polarized superfiuid. Here condensa- 
tion of Cooper pairs takes place at one or more non-zero mo- 
menta q^. These cold gases possess a remarkable flexibility 
in which they can be polarized as well as studied with vari- 
able attractive interaction. This provides an additional mecha- 
nism for (possibly) tuning in the various LOFF phases as one 
changes the s-wave two-body scattering length a from posi- 
tive in the Bose-Einstein condensation (BEC) regime to nega- 
tive (BCS). Thus far, experiments iflSEIl have focused on the 
unitary scattering regime, midway from BCS to BEC. While 
Bogoliubov-de Gennes (BdG) based theories Il9l fl0ll for a trap 
at unitarity suggest that the ground state is generalized LOFF, 
homogeneous studies full conclude that LOFF1 at least, is 
confined to a sliver near the BCS endpoint. A critical com- 
ponent which needs to be injected into this controversy is the 
nature of the stability criteria 

O Q [3 

which is, simi- 
larly, under lively debate. 

The goal of this paper is to clarify these issues by address- 
ing BCS-BEC crossover at general T. We focus on one partic- 
ular member of the LOFF class- corresponding to the single 
plane wave LOFF state, hereafter called LOFF1, in a homoge- 
neous system. We characterize the "existence regime" (where 
solutions exist) and the "stability regime" (where solutions are 
stable) in a series of phase diagrams. Unlike Ref. 11611 our 
calculations do not automatically incorporate first order tran- 
sitions (from a single phase state) to a phase-separated state. 
This phase separation is not relevant to all sub-disciplines |6J] 
and in future applications to a trap, one has also to include 
surface energy terms which are difficult to estimate. More- 
over, phase separation has not been included in the Bogoli- 
ubov deGennes (BdG) based theories flOII which we want 
to understand here. 

However, with phase separation included, our T — results 
are rather similar to those found in Ref. 1 16]. Here, we empha- 
size finite T, and we implement a highly numerical procedure 
to solve all coupled equations directly at fixed total particle 
number N a . Even in the absence of alternative phase separa- 
tion states, the stable LOFF1 state is primarily restricted to a 



regime near the BCS endpoint, although it does overlap uni- 
tarity for a narrow range of high polarizations. We show that 
the LOFF1 existence regime is considerably broader and is di- 
rectly associated with the phase space region where there 
is negative superfiuid density Il2ul8ll in the q = or Sarma 
state fisSll . Here it is likely that a LOFF phase of one form 
or another will be stable, although it may be something more 
complex OS than LOFF1. 

Our central phase diagram in the T vs p plane should be of 
particular interest to experimentalists who are currently creat- 
ing plots of this nature. It can be contrasted with that obtained 
in Ref. 112 ill in which temperature was introduced in a fash- 
ion following the original Nozieres-Schmitt-Rink (NSR) 1E2I1 
scheme. Here, unlike Ref. ll2lll we choose to include T in a 
manner which is fully consistent with the very extensive lit- 
erature 13, [l0|, InJ I12II on the ground state of these polarized 
superfluids. 

We introduce T ^ following a T-matrix scheme, and re- 
strict our attention to the superfiuid phase. This T-matrix rep- 
resents the propagator for non-condensed pairs and is given 
by t^ 1 (P) = U^ 1 + x{P) where x 1S the pair susceptibil- 
ity and U < is the pairing interaction strength. For atomic 
Fermi gases, we assume an s-wave contact interaction. There 
are strong similarities between t(P) and the Hartree-Fock ap- 
proximation to the particle-hole susceptibility which involves 
the usual Lindhard function and on-site repulsion. These 
Hartree-Fock theories are used to establish whether ferro- or 
antiferromagnetic order will arise. Here we consider a very 
similar competition between Sarma and LOFF1 states. The 
relevant x(P) necessarily involves the self consistently deter- 
mined fermionic gap parameter A(T) and chemical potential 
jtt(T). Importantly at and below T c the chemical potential for 
the pairs (Mpair) must be zero and this BEC condition on t (P), 
thereby, determines A(T). 

The pair susceptibility for LOFF1 condensates in which 
momentum k pairs with k + q, (for, as yet undetermined 
q) may readily be written down 12011 . We first introduce the 
fermionic chemical potentials /ij and fj,^ for the two spin 
states, /i = (//f + Hi)/2 and h = (/ij — /ij.)/2, and = 
k 2 /2m, £k = ek — /i, where \i a is the chemcal potential for 

spin a =t» 4- L is useful to also define i?k q = ^/?kq + 

with £ ka = (£ k + £k-q)/2 and <5e k = (ek - £k-q)/2. As 
in Ref. 1 23], we set the volume V = 1, H = ks = L and 
P = (ifti, p), where tti — 2lnT is an even Matsubara fre- 
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quency. Then the pair susceptibility x(P) at the mean field 
level, after analytical continuation iO; — * fi + «0 + , is given 
by 



/(gfcg + fek) + /(gp-k) - 1 



p-k 



/(^kq-^k)" 



,,2 /(£p-k 

k - Cp-k + (^kq - Se k ) 



(1) 



which, as q — > 0, goes over smoothly to its counterpart in 
the Sarma phase. Here the coherence factors u 2 ., = (1 ± 
&q/S kq )/2 and we define f(x) = [f(x - h) + f(x + h)]/2, 
where f(x) is the Fermi distribution function. The BEC con- 
dition, U^ 1 + x(0, q) = 0, leads to the gap equation which, 
when written in terms of the scattering length a, is of the stan- 
dard form in the literature 
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where we define f(x) = [f(x — h + Se\ ti ) + f(x + h — Se] si )]/2. 

The number equations also depend on the quantity x 
through a self energy involving t(P); we show the details 
elsewhere 1 20] but summarize our final results which yield 
the standard equations in the literature 
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(3) 
(4) 



where n — n-^ + n± and Sn = rtf — nj,, and the polarization 
p = 5n/n. Here we have defined Sf(x) = [f(x — h+ <5e k ) — 
f(x + h — 5e k )]. Finally, we determine q by imposing an 

extremal condition on the pair susceptibility, 9x gp P ^ 

0. This condition, which turns out to be equivalent to requiring 
that there be no net current in this LOFF1 state, is given by 
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When this equation has a solution at q ^ we have a LOFF1 
phase. There will always be a co-existing solution of the 
Sarma-type with q = 0. 

At T ^ 0, the parameter A(T) contains the contribution 
from both condensed (sc) and non-condensed (pg) pairs. We 
can show quite generally that below T c , A 2 (T) = A 2 C (T) + 
(T), where 



A 2 

pg 



A* (T) = Z 



_1 E 

p 



(6) 



and b(x) is the Bose distribution function. Here the pair dis- 
persion is found to be fi p w (p — q) 2 /2M*. Analytical ex- 
pressions for M* and Z are possible via an expansion of x in 
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Figure 1: LOFF1 wavevector q as a function ofl//cFfflatT = and 
for various p. Beyond the turning point, no LOFF1 state exists. The 
inset shows T c and mean field T C MF = T* (pair formation tempera- 
ture) at unitarity over respective stability regimes. 



small (p — q), 

X(fi,p)-X(0,q) 
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2M* 
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where Z = 



nil 



and 



l 

2M 



fl=0,p=q 



a=o, P=q — 6Z 9 p 2 

The quantity x contains everything one needs to know about 
zero as well as finite T. And our results for T = reduce to 
the standard equations in the literature. 

To demonstrate that a given LOFF1 solution to our self con- 
sistent equations is stable, we introduce an effective thermo- 
dynamic potential for this superfluid state 

^ = -^+£{(£kq-£kq) 
k 

- Tln(l + exp[-(£ kq - h + Se k )/T}) 

- Tln(l + cxp[-(£ kq + h - 5e k )/T})} . (8) 

It is straightforward to verify that the above gap, number and 
zero-current equations are consistent with the variational con- 
ditions 

— -Sn- — 
dh ' dq 

The stability condition requires that the symmetric number 
susceptibility matrix 
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be positive definite fiaHai.Here^ = ^+f£4+§f.JL, 
with x = n, h. To evaluate this matrix we note 
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Figure 2: (Color online) Phase diagram in the p vs 1/kpa plane for T/Tf = 0, 0.05, and 0.1. The dotted region shows where LOFF1 state is 
unstable but some form of LOFF phase may be stable. The (yellow) light shaded region indicates the stable LOFF1 superfluid and the (cyan) 
darker shaded region the stable LOFF1 normal state. 



where dA/d/J,, dq/dfi, dA/dh, and dq/dh can be easily de- 
rived by differentiating Eqs. (0 and (0 with respect to /j and 
h. It can be shown that the positive definiteness of M is equiv- 
alent to 

mw-(m^) >0 ' (10) 

Equations (|2j-(|6} were solved numerically to establish the 
existence regime for LOFF1 solutions. Figure^shows a plot 
of the behavior of the LOFF1 wave vector q as a function 
of 1/fcpa, and at T = 0, for various polarizations p. Here 
Ef = UbTf = ti 2 k F /2m is defined as the Fermi energy of 
an unpolarized, noninteracting Fermi gas of density n. For 
each value of p a turning point, {l/kFa) max , is visible be- 
yond which we can not find LOFF1 solutions. At finite T the 
analogous curves (not shown) quickly become monotonically 
decreasing so that (l/fci?a) maa; corresponds to q = 0; thus 
the system smoothly transforms to the Sarma state. 

The inset in Fig.n snows the behavior of T c and its mean 
field counterpart T* = T™ F as a function of 1/fcpa near 
unitarity and for p = 0.8. It can be seen that there is a con- 
siderable difference between T c and T^ IF showing that pair 
fluctuation effects (via A pg ^ 0) are very important in the 
LOFF1 phase, just as seen elsewhere 11811 . In contrast to the 
behavior of the Sarma phase at unitarity, in the LOFF1 state, 
superfluidity extends over the range of temperatures from T c 
down to T = 0. 

We turn to Fig. [2] beginning with the left panel (T = 
0) and the nearly vertical line which is determined from 
(l/kFa)max- This provides a bound on the existence region 
for our numerically obtained LOFF1 solutions. Essentially on 
top of this nearly vertical line is the locus of points to the left 
of which the superfluid density for the Sarma state nf ar " la (0) 
is negative. That these two lines coincide reinforces earlier 
work 1 17]: providing one considers a second order transition 
between generalized LOFF and Sarma states, the boundary 
line for the existence of all LOFF states is determined by 
n Sarma _ q m( ji C ate by the dotted background in Fig. [2] 
where we have the possibility of superfluidity with multiple 
nonzero q's. It should be noted that this existence regime for 



generalized LOFF states is relatively wide at T = 0, impor- 
tantly including unitarity. The shaded region in Fig.|2]results 
from applying the stability criteria associated with the posi- 
tivity of the matrix M in Eq. (|9jl. The non-vertical line to the 
right in this first panel marks the onset point for a stable Sarma 
phase. 

This phase diagram evolves with temperature as shown by 
the other two panels in Fig. [2] In the middle and right panels 
we have distinguished between normal and superfluid LOFF1 
states by using darker and lighter shaded regions respectively. 
At the highest T, for the panel on the right, there is only a 
sliver of stable LOFF1 which corresponds to a normal (pseu- 
dogapped) phase. It should be noted that the size of the ex- 
istence region for generalized LOFF solutions quickly de- 
creases as temperature is raised. This is related to the fact 
that the nf arma rapidly becomes non-negative as T increases 
from zero. 

Figure [3] represents a particularly convenient way of pre- 
senting our results. For the trapped case, experimental studies 
01 HI 01 are in the process of mapping out this phase diagram 
in the p-T plane at unitarity. From left to right, the three 
panels correspond to unitarity, and the BCS and BEC sides 
(close to resonance), respectively. Also shown here is the re- 
gion where we have a stable Sarma state. This appears only at 
intermediate temperatures |18], when the superfluid density 
(which is negative at T = for all 3 cases) is driven posi- 
tive. Note that our p-T phase diagram is different from that in 
Ref . 1 2 1 ] , which is based on a different but unspecifed ground 
state. Indicated in all three panels are the (dotted) regions 
where generalized LOFF states may exist, and the (shaded) 
regimes where the LOFF1 phase is stable. As in Fig. [2] 
the light shaded region corresponds to the superfluid LOFF1 
phase, and the dark shaded region to the normal LOFF1 phase 
with a pseudogap. It is evident that, for these 1/kpa values, 
a stable LOFF1 phase exists only at relatively high p and low 
T. This is to be contrasted with the stable Sarma superfluid 
which exists only at low p and intermediate T. Using this 
figure, one can compare the transition temperatures for the 
LOFF1 and Sarma phases. For the latter, T c is read off as the 
upper transition temperature in the 3 plots. (The lower T c is 
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Figure 3: (Color online) Phase diagram in the p — T plane for unitary (left), near BCS (middle) and near-BEC (right). The dotted region shows 
where LOFF1 state is unstable but stable generalized LOFF phases may in principle exist. The (yellow) light shaded region is the LOFF1 
superfluid, and the (cyan) darker shaded region is the LOFF1 normal state. 



where n s 8 arma (T) chan ges sign and the order parameter van- 
ishes). It can be seen that the transition temperatures for the 
LOFF1 phase are very low compared to their counterparts in 
the Sarma phase. 

In summary, in this paper we have addressed homogeneous 
systems and mapped out the LOFF1 phase diagram at general 
T. We have shown that LOFF1 phases are more stable near 
unitarity at sufficiently high p than alternative Sarma states. 
Since it is generally expected |7] that LOFF states with mul- 
tiple values of q are more stable than the simplest LOFF1 
states, we argue that quite possibly there exist stable general- 
ized LOFF ground states throughout most (but perhaps not all) 
of the dotted regions shown in Figs. 2 and 3; this corresponds 
to where the Sarma phase has negative superfluid density. Our 
results support previous BdG-based approaches |9, 10] which 



argue that the ground state at unitarity is a generalized LOFF 
phase. In a trap configuration, one might expect that, since we 
find the LOFF1 phase is stable at relatively high p, generalized 
LOFF phases should appear in the neighborhood of the con- 
densate edge, as also found in earlier work |9, 10]. However, 
in contrast to Ref. 0], we have addressed systematic LOFF1 
stability criteria and, moreover, find that the size of the sta- 
bility region for the LOFF1 state and the size of the existence 
region for general LOFF phases very quickly diminish with 
T. Finally, our self consistent calculations indicate very low 
T c in the LOFF1 state as compared with previous estimates in 
the literature |9], which ignored the effects of non-condensed 
pairs. 
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